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We study boron, carbon, nitrogen and oxygen isotopes with a newly constructed shell-model 
Hamiltonian developed from monopole-based-universal interaction (Vmu)- The present Hamiltonian 
can reproduce well the ground-state energies, energy levels, electric quadrupole properties and spin 
properties of these nuclei in full psd model space including (0 — 3)ftu excitations. Especially, it 
correctly describes the drip lines of carbon and oxygen isotopes and the spins of the ground states 
of 1C, B and 18 N while some former interactions such as WBP and WBT fail. We point out that the 
inclusion of 2hw excitations is important in reproducing some of these properties. In the present 
(0 + 2)hw calculations small but constant E2 effective charges appear to work quite well. As the 
inclusion of the 2Tiw model space makes rather minor change, this seems to be related to the smallness 
of 4 He core. Similarly, the spin g factors are very close to free values. The applicability of tensor 
and spin-orbit forces in free space, which are taken in the present Hamiltonian, is examined in shell 
model calculations. 



PACS numbers: 21.60.Cs, 21.10.-k, 23.20.-g, 23.40.-s 

I. INTRODUCTION 

The existence of the unexpected doubly magic nucleus 
24 shows the exotic property of drip-line nuclei, that is, 
the change of magic numbers far from the stability [l|. 
One of the aims of theoretical works on nuclear structure 
is to describe both stable nuclei and nuclei far from the 
stability in a unified framework. In shell-model studies, 
for many of existing conventional interactions, it is diffi- 
cult to reproduce simultaneously the drip lines of carbon 
and oxygen isotopes as well as some other properties such 
as energies of 2^~ states and B(E2). From a microscopic 
study, the inclusion of effects of three-body force is im- 
portant in describing the drip line of oxygen isotopes Q ■ 
It is emergent to construct new shell-model interactions 
applicable from the /3-stability line to the drip lines. 

The realistic nuclcon-nuclcon (NN) interactions need 
to be renormalizcd when applying to shell-model cal- 
culations because of the short-range correlation and in- 
medium effect NN interaction is composed of three 
components, central force, spin-orbit force and tensor 
force. Recent studies show that the monopole compo- 
nents of tensor force barely change after the renormal- 
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ization and that the multipole components also change 
little [H H[. Based on these studies, a monopole-based- 
universal interaction (Vmu) including the bare ir + p ten- 
sor force is introduced to describe the shell evolution Q ■ 
As this Vmu is constructed based on monopole proper- 
ties, it requires examination as to whether Vmu can be 
used in actual shell-model calculations or not. In this pa- 
per we try to apply the Vmu to shell-model calculations 
in psd model space. 

In the psd region, several effective interactions have 
been introduced in shell-model calculations, such as PS- 
DMK @, WBT 0, WBP Q and SFO [|. PSDMK, 
WBT and WBP interactions are all constructed in (0 — 
l)hu> model space, which means that — 1 nucleons are 
allowed to be excited from p shell to sd shell. The mix- 
ing between (0 — l)hu! states and (2 — 3)hu; states is not 
considered in the fitting of the interaction. SFO, which 
includes the (2 — 3)Huj states, concentrates mostly on the 
spin properties such as magnetic moments and Gamow- 
Teller transitions. Up to now, the (pp\V\sdsd) matrix ele- 
ments, which represent the interaction between (0 — l)huj 
states and (2 — 3)fiw states, have not been well studied. 
In Ref. @, the tensor part of the (psd\V\psd) matrix 
elements is taken to be that of the it + p meson exchange 
potential and spin properties of C isotopes have been 
studied. Recently, the study of microscopic derivation of 
the effective interaction for the shell model in two major 
shells is in progress [IFJ . It would be interesting to apply 
them to shell-model calculations in future. 

In this paper we try to construct the effective in- 
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teraction in the psd space based on Vmu to describe 
ground-state energies, energy levels, electric quadrupole 
properties and spin properties. The {psd\V\psd) and 
(pp\ V\sdsd) matrix elements are obtained based on Vmu 
while phenomenological effective interactions are used 
for the p-shell and sd-shell parts to maintain the good 
description of the phenomenology by these interactions. 
Microscopic interactions have been obtained based on G- 
matrix method with medium modification fill ] , similarity 
renormalization group (SRG) method [l2j and coupled- 
cluster method [13| ]. While they produce interesting re- 
sults, fully microscopic calculations have not been suc- 
cessful, as far as a good agreement to experiment is con- 
cerned. We restrict here to a more phenomenological 
approach based on Vmu to study the spectroscopic prop- 
erties of the nuclei to be discussed. 

In the next section, we introduce a new Hamiltonian. 
Coulomb correction and center-of-mass correction are 
discussed in Sec. IIIII and Sec. IIV1 respectively. In Sec.lVl 
we discuss the ground-state energies and energy levels. 
We present the results of electric quadrupole properties 
and spin properties in Sec. IVII and Sec. IVII1 respectively. 
A summary is given in Sec. IVIIII 



II. HAMILTONIAN 

The present Hamiltonian is developed from Vmu, 
SFO and SDPF-M Q. The two-body matrix ele- 
ments (TBME) are constructed as follows, (pp\V\pp) 
from SFO, (sdsd\V\sdsd) from SDPF-M, (psd\V\psd) 
and (pp\V\sdsd) from Vmu phis spin-orbit force. In 
the (pp\V\pp) matrix elements, we reduce the strength 
of the monopole term (pi/2P3/2|^|pi/2P3/2)t=o by 0.5 
MeV from SFO. This will improve the description of the 
ground-state energies of these nuclei. The (sdsd\V \sdsd) 
matrix elements in the present Hamiltonian are the same 
as SDPF-M. In earlier interactions, such as WBP and 
WBT, the matrix elements (pp\V\sdsd) are not consid- 
ered in the fitting procedure. The strength of the in- 
teraction in (pp\V\sdsd) in WBP and WBT is the same 
as in (psd\V\psd) in WBP. In the present interaction, 
strengths of these two parts of the interaction are not 
taken to be the same. The Vmu includes Gaussian type 
central force and it + p tensor force. We use M3Y 15 1 
force for the spin-orbit force. We keep the spin-orbit and 
tensor forces unchanged. The form of the interactions in 
the matrix elements of (psd\ V\psd) and {pp\V\sdsd) is as 
follows, 



V Vcentral V S pi n — orbit 

(M3Y) + V t 

elisor fr + p), (1) 

with V cen trai being, 

Vcentral = ^ fs ,T P S,TeXp(- (r / fl) 2 ) , (2) 
S,T 

where S (T) means spin (isospin), Ps.t is the projec- 
tion operator on S, T channel, r and /x are distance 



between two nucleons and Gaussian parameter, respec- 
tively. fsT is the strength of the central force. In the 
original V M u, fo,o = fi.o = -166 MeV, / ,i = 0.6/ ,o 
and /i i = 0.8/o,o 3- In the present study, we reduce 
the central force in (psd\V\psd) and (pp\V\sdsd) matrix 
elements by factors 0.85 and 0.55 from the original Vmu, 
respectively The final interaction in the (psd\V\psd) 
((pp\V\sdsd)) matrix elements is 



V = 0.85(0.55)F ce „ tro ; + 

* spin — orbit 

(M3Y) + V t en sor 

(tt + p). 



(3) 



Notice that the spin-orbit force and the tensor force are 
kept unchanged. The TBME are calculated with har- 
monic oscillator parameter fku = 45 A" 1 / 3 — 25A -2 / 3 
where A = 18 which is the average mass number of the 
investigated nuclei from 10 B to 26 0. The sd-shcll single- 
particle energies (SPE's) in SDPF-M are e drj/2 = -3.95 
MeV, e d3/2 = 1.65 MeV and e Sl/2 = -3.16 MeV, which 
takes 16 O as the core (l4| . In the present shell-model 
calculations, 4 He is chosen as the core, thus the sd-shell 
SPE's in the present Hamiltonian should be adjusted to 
give the same one-particle excitation energies of 17 as 
in SDPF-M. The adjusted SPE's are e d5/2 = 8.01 MeV, 
e d3/2 = 10.11 MeV and e Sl/2 = 2.11 MeV. The p-shell 
SPE's are obtained based on SFO but with slight changes 
by fitting the ground-state energies of the studied nu- 
clei and related levels such as the 1/2]" state in n B and 
3/2j7 state in 13 C. We obtain e p ,„ = 1.05 MeV and 
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P3/2 

= 5.30 MeV. The detailed TBME of the present 
Hamiltonian can be obtained by contacting the authors. 

We compare the TBME of the present Hamiltonian 
with those of WBT and WBP in Fig. [TJ The TBME of 
central, spin-orbit and tensor interactions are also pre- 
sented by the spin-tensor decomposition method [16j |. 
The sd and ppsdsd parts of WBT and WBP are the same 
between the two. So we show only the WBT result in 
these two parts. The sd part of the present interaction 
is from SDPF-M which is modified from USD (the same 
as the sd part of WBT) interaction. There is not much 
difference between WBT and the present interaction in 
the sd part. In the p part, all these three interactions, 
present, WBT and WBP, are fitted to low-lying levels of 
the p-shell nuclei. The difference among these three in- 
teractions is not large except for the tensor force. In the 
psd and ppsdsd parts of the interaction, the deviation of 
the present interaction from WBT (WBP) turns out to 
be larger. The central force of the present interaction in 
the psd part is 0.85Vmu- We find that this strength is 
proper as the number of points above the diagonal line is 
close to that below the line as shown in the Fig.QJc). It is 
interesting that the spin-orbit interaction of the present 
interaction is very similar to that of WBP interaction in 
both the psd and ppsdsd parts. In these two parts of 
the interaction, WBP has 10 parameters for the poten- 
tial fitting while the present interaction is taken from the 
M3Y potential. Quite similar results between WBP and 
the present interaction indicate that the spin-orbit force 
is rather well determined compared to the central force. 
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FIG. 1: (Color online) Total TBME, TBME of central force, spin-orbit force and tensor force in each part of interaction: (a) 
(sdsd\V\sdsd) , (b) (pp\V\pp), (c) (psd\V\psd) and (d) (pp\V\sdsd) . 



Figure [5] presents the monopole terms of the interac- 
tions and their spin-tensor components. Monopole term 
is a weighted average of TBME for orbits j and j' [13, EH , 

v t = EA2J + i)(jf\v\jf)j,T m 
E,/(2J + i) ' [) 



The monopole terms are presented in three groups, pp, 
sdsd and psd in each picture. In each group the cen- 
tral monopole is relatively flat compared with the total 
monopole. The total interaction can be recognized as a 
global central force plus other staggers. The T = cen- 
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FIG. 2: (Color online) Monopole terms of total, central, spin-orbit and tensor forces in each isospin channel: (a) T = and (b) 
T = 1. 
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tral monopole is the most attractive among all these six 
central, spin-orbit and tensor monopoles. The nuclear 
binding energy comes mostly from this interaction. Both 
T = and T = 1 spin-orbit monopoles of the present in- 
teraction are very close to those of WBP. This is consis- 
tent with the analysis of the spin-orbit part of the TBME. 
Comparing with WBT, the present spin-orbit monopoles 
are also not much different. The present tensor force is 
stronger than WBT and WBP in T = channel, more 
attractive in (j>i/2^5/2|^bi/2^5/2) an d more repulsive in 
(^3/2^5/21^1^3/2^5/2)- In the sd region of nuclei, this ef- 
fect of the tensor force is canceled as the pi/ 2 and p 3 / 2 
orbits are fully occupied [lj|. Going to psd region, such 
as neutron rich boron, carbon and nitrogen isotopes, the 
opposite sign of the monopoles of the tensor force turns 
out to be important. 



III. COULOMB CORRECTION 

In the shell-model study, Coulomb interaction is not 
included in many cases, in order to keep the isospin sym- 
metry. When we compare the ground-state energies be- 
tween theoretical results and observed values, Coulomb 
correction is needed. Present calculations in psd model 
space do not include the ground-state energies of 4 He, 
E( 4 He), which also needs to be removed. The total cor- 
rection is as follow, 



E, 



correction 



Eco 



E( 4 He), (5) 



where E Coulomb and E correction are the energy of 
Coulomb correction and the ground-state energy after 
the correction, respectively. E( 4 He) = —28.296 MeV. 
Ecouiomb is calculated through similar method used in 
the construction of WBT and WBP interactions 0- 
We calculate the energy difference of mirror nuclei near 
N = Z where the observed ground-state energies are 
taken from [2Cj ■ This Ecouiomb is dependent only on Z in 
our calculation. Ecouiomb = 1-075 [Z = 3), 2.720 (Z 
4 . 4.593 (Z = 5), 7.368 (Z = 6), 10.248 (Z = 
7), 13.854 [Z = 8) MeV. 
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IV. CENTER-OF-MASS CORRECTION 

As our calculation is done in two major shells, we need 
center-of-mass (cm.) correction to remove the spurious 
components which come from the cm. motion. We use 
the method suggested by Gloeckner and Lawson [2l[ . In 
the calculations, Hamiltonian is H' = Hsm + /3J? C . TO ., 
where Hsm and H c . m . are original and cm. Hamiltoni- 
ans, respectively. If (3 is large enough, the effect of the 
cm. motion is small enough in low lying states. Fig- 
ure [3] indicates some physical quantities of 16 C to check 
whether this method works or not, and how large /3 is 
needed. 

We find that the number of nuclcons in the sii-shell, 
B(E2;0+ -> 2+) or the energy of 2f in 16 C hardly change 



FIG. 3: (Color online) Effects of the center-of-mass correc- 
tions in four physical quantities: (a) nucleon number in sd- 
shell, (b) B(£'2;0+ 2+) value, (c) E(2+) and (d) ground- 
state energy of 16 C. Their values change as the function of 
P. 



when /3 changes. The ground-state energy of 16 C changes 
quickly when (3 is small. For j3 > 10, it becomes almost 
flat. We use j3 = 10 in the following calculations. 
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FIG. 4: (Color online) Ground-state energies of boron, car- 
bon, nitrogen and oxygen isotopes. Experimental values are 
taken from Ref. [2y| 
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FIG. 5: (Color online) Two neutron separation energies, S2 n , 
of boron, carbon, nitrogen and oxygen isotopes. Experimental 
values are taken from Ref. [201 



V. GROUND-STATE ENERGY AND ENERGY 
LEVEL 

The nuclei 22 C, 23 N and 24 are the last bound nuclei 
in neutron rich side of C, N and O isotopes 



22(. The 



neutron-drip lines in elements beyond oxygen is not de- 
termined yet [H]. In WBT and WBP, 22 C is unbound 
and 26 O is bound. The present Hamiltonian improves the 
description of drip lines of C and O isotopes. Figures @] 
to [6] present the ground-state energies as well as one and 
two neutron separation energies S n and 52„, for B, C, 
N and O isotopes. In O isotopes, the WBT and WBP 
have the same result as their sd parts. From S n and 




10 - 



> 

CD 



if) 



-10 



N 



(b) 



■ Expt V \/\ 



WBT 

WBP 

present 



10 



15 



A 



20 



25 



FIG. 6: (Color online) One neutron separation energies, S n of 
boron, carbon, nitrogen and oxygen isotopes. Experimental 
values are taken from Ref. I20II. 
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FIG. 7: (Color online) Comparison of S„ and S2n under three 
different values of the center-of-mass parameter, ,8 = 1, 10 
and 20. 
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FIG. 8: Energy levels for boron and light carbon isotopes, obtained in the present, SFO and WBP calculations, compared with 
experimental data [25j - l3(il |. 
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FIG. 9: Similar to Fig. but for heavier carbon and light nitrogen isotopes. 
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FIG. 10: Similar to Fig. [8] but for heavier nitrogen and light oxygen isotopes. 
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FIG. 11: Similar to Fig. El but for heavier oxygen isotopes. 
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S2n of O isotope, one can see that both WBT and the 
present Hamiltonian predict 25 O to be unbound consis- 
tently with experiment H3- The 26 is about 1.2 MeV 
unbound in the present result and 1.0 MeV bound in 
WBT. The positive S n value in 26 indicates that 26 is 
one-neutron bound but two-neutron unbound. In N iso- 
topes, all these three interactions can reproduce that 23 N 
is bound and 24 N and 25 N are unbound. 21 C is unbound 
experimentally [22| and also unbound in the calculations 
by all these three interactions. 22 C is 0.1 MeV bound in 
the present result and 0.2 MeV and 0.6 MeV unbound 
in WBT and WBP, respectively. 16 B is 40(60) keV un- 
bound [H. It is 144 keV bound in WBP and 65 keV and 
153 keV unbound in WBT and the present Hamiltonian, 
respectively. All Hamiltonians succeed in describing un- 
bound 18 B. In the experiment [22|, WBP and present 
Hamiltonian, it is one-neutron unbound. But in WBT 
it is both one and two-neutron unbound. 19 B, which is 
experimentally bound [23], is unbound with 160 keV, 381 
keV and 538 keV in the present Hamiltonian, WBT and 
WBP, respectively. 

Here we briefly summarize the descriptions of drip lines 
by these three Hamiltonians. The present Hamiltonian is 
successful in describing all drip-line nuclei except for 19 B 
while WBT fails in 26 , 22 C, 18 B and 19 B and WBP fails 
in 26 0, 22 C, 16 B and 19 B. One reason that the present 
interaction improves the description of drip lines is the 
inclusion of the mixing between OHuj and 2hw configura- 
tions. WBT and WBP have mass-dependent term in the 
sd-shell 0- Going from 18 O to 28 O, the sd shell interac- 
tion decreases, which makes the nuclei less binding. We 
find that the mixing between Ohu and 2hjj states has a 
similar effect. Partial effect of mass dependence there- 
fore comes from the mixing between OHlj and 2Huj states 
which is not included in WBT and WBP. Wc will discuss 
more about the contribution of 2Hlo states later. 

In order to see if the prediction on the neutron drip line 
is sensitive to the center-of-mass parameter, (3, we have 
made calculations with assuming three different values 
of the parameter. Figure [7] displays the calculations of 
one- and two-neutron separation energies under differ- 
ent (3 values, showing that the value of /3 = 10 used in 
the present work is large enough to remove the spuri- 
ous center-of-mass components. For example, the neu- 
tron separation energies of neutron-rich oxygen isotopes 
change about 100 keV when increasing (3 from 1 to 10, 
while the separation-energy variation is about 20 keV 
when increasing j3 from 10 to 20. It is consistent to dis- 
cussions in Sec. IIVI that the physical properties are well 
convergent when /3 = 10. 

Figures [81 to [Til present the energy levels of B, C, N, 
and O isotopes. The agreement between experiment and 
the present work is fairly good. Especially for 10 B and 
18 N, we can reproduce the spins of the ground states of 
these two nuclei while WBP and WBT fail. WBT also 
fails in describing the spins of the ground states of 17 C, 
19 C and 16 N. Wc only show WBP results here as WBP 
is similar to WBT and is better in the description of the 
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spins of the ground states. The ground states of 16 N 
and 22 N are about 100 keV higher in the present interac- 
tion. The SFO can also reproduce the spin of the ground 
state of 10 B and nuclei nearby. But it fails in some neu- 
tron rich nuclei such as 19 O and 21 0. This is because 
the (sdsd\V\sdsd) part of SFO is from renormalized G- 
matrix Q. The interaction is too attractive without the 
contribution of three-body forces 0. The first l/2 + and 
5/2+ states in 19 0, 21 and 23 indicate that the neu- 
tron lsx/2 orbit is too low compared with neutron 0d 5 / 2 
orbit in SFO. This situation can be improved by using 
effective interactions such as SDPF-M or including the 
contribution of three-body forces @], for instance. 

The energy difference between first 3 + and 1 + in 10 B 
can be reproduced well by both SFO and the present 
Hamiltonian. This is partly because the Opi/2 orbit is 
much higher than Op^/2 orbit and partly because the 
strength of (pp\V\sdsd) is chosen properly which will 
be discussed later. For 10 B and 11 B, an ab initio no- 
core shell model calculation based on chiral perturba- 
tion theory showed that the inclusion of the three-body 
force is necessary to reproduce the ground-state spins 
[23l ]. In the present Hamiltonian, the phenomenologi- 
cal effective two-body interaction is mostly obtained by 
fitting experimental data. Therefore, the effective inter- 
action obtained thus includes, at least partly, the three- 
body effect. Our calculations show this equivalence. In 
Ref. , it was pointed out that ab initio interaction with- 
out three-body force cannot reproduce the neutron drip 
line of oxygen isotopes, while ab initio interaction with 
three-body force or phenomenological two-body interac- 
tions may describe the drip line. 

The v(sd) 3 configuration shows different structure in 
N = 11 isotones from 17 C to 19 0. i>(0d 5/2 ) 3 can 
couple to J = 5/2 with seniority v = 1 or couple 
to J = 3/2 with seniority v = 3. The structure of 
is(sd) 3 as well as the low lying states in iV = 11 iso- 
tones is a subtle problem because of these two configura- 
tions together with v{0d 5 / 2 ) 2 (ls 1 / 2 ) 1 , v(0d 5/2 ) 1 (ls 1 / 2 ) 2 
and many small components. For the first time, the 
present Hamiltonian in full psd model space reproduces 
the low lying states in all these three nuclei, 17 C, 18 N 
and 19 0. {{0d 5 /2) 2 \V\(0d 5 / 2 ) 2 ) j=o t=i paring interaction 
contributes to this good agreement because this paring 
is reduced in SDPF-M which will make the v = 1 state 
less bound and keep v — 3 state unchanged. As one can 
see, 2+ state in 18 N, 5/2+ states in 19 and 17 C in the 
present results become higher compared with those in 
WBP results. Other matrix elements also contribute to 
this subtle problem. We will discuss the contribution of 
(pp\V\sdsd) in 17 C later. 

The WBP and WBT results show more expanded en- 
ergy levels compared with observed energy levels in C 
and N isotopes [2^, [2fj| . This can be improved by reduc- 
ing neutron-neutron interactions by 25% (for C isotopes) 
or 12.5% (for N isotopes) in the sd-shell in WBP and 
WBT [25], |2(| . The spectra of the present interaction are 
not so expanded as in WBP and WBT for C, N and O 



isotopes. 

In the present work, energy levels of unnatural parity 
state are not fully considered. One reason is that exper- 
imental data of these energy levels are not much avail- 
able in neutron rich nuclei. Another reason is that the 
dimension of the calculation increases quickly when in- 
cluding Shio components. We can improve the description 
of these unnatural parity states with more experimental 
data and more advanced computers in the future. 

The strength of the interaction in the (pp\ V\sdsd) ma- 
trix elements is not determined in PSDMK, WBP and 
WBT and not fully considered in SFO. In the present in- 
teraction, the strength of the central part of (pp\V\sdsd) 
is 55% of Vmu- We will show some examples that the 
(pp\V\sdsd) matrix elements are important in describ- 
ing the nuclei being studied. The total wave function of 
a nucleus can be written as $ = a^>(0hu) + b^(2hu). 
Figure [T2l shows the probability, b 2 , of 2hw component. 
It is clear that the probability b 2 is very sensitive to 
neutron numbers. When the neutron number increases 
from 8 to 15, the value of b 2 decreases except a singu- 
lar point 17 C. In WBP and WBT, the sd part includes 
the mass-dependent term (18/A) - 3 0). Only with this 
effect, WBP and WBT can reproduce well the ground- 
state energies of these nuclei. The mass-dependent term 
is needed for calculations of nuclei in a large mass range 
because the nuclear force is related to the radii of nuclei 
as well as the nucleon number A. But in a range of nuclei 
with small mass numbers, the effect of mass dependence 
is not obvious when we include 2hio components. In the 
present Hamiltonian, we can reproduce well ground-state 
energies, separation energies and energy levels of B, C, N 
and O isotopes without mass dependent term. One can 
see from Fig.[l2]that the inclusion of 2fruj components will 
automatically contain a part of mass-dependent effects. 
More works are needed to study the mass dependent ef- 
fects in light nuclei. 

The (pp\V\sdsd) matrix elements are also important 
for energy levels in certain nuclei. Figure [T3] shows the 
dependence of the energy levels in 10 B and 17 C on the in- 
teraction. Energy differences, such as difference between 
3+ and if in 10 B and that between 3/2+ and 5/2+ in 17 C, 
are very sensitive to the strength of {pp\V\sdsd}. The en- 
ergy difference between 3 + and + in 10 B, on the other 
hand, is hardly changed when the strength of the cen- 
tral part of (pp\V\sdsd) is changed by 60% of Vmu- The 
above observations suggest that the contribution of 2fku 
components is not only A dependent but also state de- 
pendent. The 2Hu> components are 4.3%, 16.0% and 6.0% 
in 3 + , 1 + and + states in 10 B, respectively. It is inter- 
esting to do systematic investigation on how (pp\V\sdsd) 
as well as 2hio or more fko components affect the ener- 
gies and effective operators, such as effective charges and 
spin g factors which will be mentioned in the next two 
sections. 
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VI. ELECTRIC QUADRUPOLE PROPERTIES 

The present Hamiltonian has been shown to be able to 
describe the energies of the psd-sheW nuclei quite well. It 
is necessary to investigate whether this interaction gives 
appropriate wave functions as well. In this section, we 
discuss the electric quadrupole properties with the use 
of the present Hamiltonian and WBP. In shell model, 
effective charges are needed because of the polarization of 
the core which is not included in the model space 
One set of effective charges, e n = 1.3 and e„ = 0.5, is 
suitable for sd-shell nuclei [31(, which means that both 
valence protons and neutrons are excited in the s<i-shell. 
For valence protons and/or neutrons locate in p shell in 
neutron-rich nuclei, this set of effective charges becomes 
invalid HHIt 

Figure [14] shows the quadrupole moments in B, C and 
N isotopes and B(E2) in Be and C isotopes with two 
sets of effective charges, one is Z, N dependent [IH and 
the other is independent of Z and N. Experimental val- 
ues are taken from Refs. (33 - l4l| . For the Z and N in- 
dependent effective charges, we obtain them by fitting 
to quadrupole moments of these nuclei except for 18 N 
and 10 B. Quadrupole moment of 18 N is not exactly de- 
termined as there are two experimental values [341 ] . In 
case of 10 B, Z- or N-indcpendcnt effective charges cannot 
describe well its quadrupole moment, as will be discussed 
later. 

The Z and N independent effective charges obtained 
for the present Hamiltonian and WBP are e p = 1.26, 
e„ = 0.21 and e p = 1.27, e„ = 0.23, respectively. We 
also get the effective charges for the present Hamiltonian 
in Oftu model space, e p = 1.25 and e„ = 0.25. The inclu- 
sion of the 2hu! model space reduces the effective charges 
a little. Both of them underestimate the quadrupole mo- 
ments in stable nuclei such as 10 B, 11 C and 12 N and 
overestimate those of the nuclei somewhat far from the 
stability-line such as 15 B and 17 B. This probably means 
that stable nuclei have stronger core polarization while 
nuclei far from the stability-line have weaker core polar- 
ization. In nuclei far from the stability-line, some valence 
nucleons are weakly bound, which will make the radial 
wave function extended farther than the well bound nu- 
cleons. The extended wave function will reduce the in- 
teraction between valence nucleons and the core. In case 
of 12 C, the results for B(E2) values with fixed effective 
charges are better than those with Z and N dependent 
effective charges. We emphasize that the smaller but con- 
stant effective charges can reproduce experimental data 
rather well in Fig. 1141 The smallncss may be explained 
as a consequence of the small core of 4 He in the present 
work. More studies on effective charges are of great in- 
terest. 

Although none of combinations of WBP or present 
Hamiltonians with either set of effective charges works 
well in the quadrupole moment of 14 B, the present Hamil- 
tonian improves the result of 14 B compared with WBP. 
We also calculate this quadrupole moment in Ohio model 



space with the present interaction. The result becomes 
worse than that we show in Fig. Q3] which is obtained in 
2hio model space. The ground state of 14 B includes 17% 
of 2huj configurations. Including more hio excitations 
may improve the result of 14 B. The 2hio configurations 
also improve the B(E2;0+ 2+) of 12 Be. The 2+ of 
12 Bc is almost a pure 2hw state, that is, with 93% of the 
2hw components in the present Hamiltonian. The Of and 
Of of 12 Be have 64% and 54% of the 2ftw components, 
respectively. Therefore, although B(E2; Of —> 2+) val- 
ues are very close to each other in WBP and the present 
results, they are contributed by different configurations 
in each calculation. In the WBP result, B(E2; Of ->■ 2|) 
of 12 Be is all from the contributions by the transition 
between p shell nucleons, especially the transition in- 
side Ops/2 proton orbit. In the present result, besides 
p shell protons, p and sd shell neutrons contribute a lot 
to B(E2; Of -)■ 2+) in 12 Be. In 12 Be, the pure p shell 
proton is not enough to reproduce the B(E2;0f — > 2+) 
value as we see in Fig. Q31 

We also try the conventional effective charges for sd 
shell, e p = 1.3 and e„ = 0.5, to calculate quadrupole 
moments and B{E2) values with the present Hamiltonian 
in Fig. [TS] It is seen clearly that this set of effective 
charges is also invalid for this new Hamiltonian. Almost 
all values are much overestimated with this set of effective 
charge. 



VII. SPIN PROPERTIES 

If two protons (neutrons) couple to a pair of angular 
momentum zero, their total magnetic moment (m.m.) is 
zero. The m.m. reflects the motion of unpaired protons 
and/or neutrons. Figure [16] presents the m.m with WBP 
and the present Hamiltonian in both Ohio and 2fku model 
spaces with 5g$ v = ±0.1 and gi e * jg s = 0.95,0.92 
and 0.90 for the present 2fiw, the present Ohio and WBP, 
respectively. 5g^ comes from the meson exchange pro- 
cesses [42|, |43| and Sg^/g^ is obtained from the \- 
square fitting of the calculated values to the experimen- 
tal ones in these nuclei. All theoretical results reproduce 
well the observed values except for a few nuclei. The 
largest deviation between calculations and experimental 
results is found in 18 0, 16 N, 14 C and 14 B. Notice that 
their E2 properties are also not well described. We do 
not show the result for 18 O in the present work. Its cal- 
culated B(E2) value is much smaller than the observed 
one. These nuclei may demand larger model space with 
Ahio or more excitations. 

In Fig. I17[ we also show the root-mean-square (r.m.s.) 
deviation between calculations and observed values, 
\m.m.(Exp.) —m.m.(Cal.)\, versus g^f^^ / g s . In a region 

5gf^^ I g s = ±0.03, the r.m.s. deviation of each result 
is very flat. Outside this region, the r.m.s. deviation in- 
creases. The minimal point for the r.m.s. deviation is 
located at gi eff) /g s = 0.95, 0.92 and 0.90 for the present 
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to be e p — 1.26, e n = 0.21 and e p = 1.27, e n = 0.23 for the present and WBP, respectively, (a) Electric quadrupole moments 
calculated with Z, N dependent effective charges, (b) B(E2) values calculated with Z, N dependent effective charges, (c) 
Electric quadrupole moments calculated with fixed effective charges, (d) B(E2) values calculated with fixed effective charges. 
All quadrupole moments are for the ground states except for 2^ in 12 C. All B(E2) values are from 0^ to 2^ except for the 
second B(E2) value in 12 Be, which is from to 2+. 



TABLE I: B(GT) values of experiment, WBP, present in both Oftu) and 2ftu) results. Experimental values are taken from the 
Ref. [H and the related references in this paper. 



2tkJ, the present Ohui and WBP, respectively. As we ex- 
pect, the quenching is weaker when we enlarge the model 
space. The quenching of the present 2fku result is rather 
weak and we may safely use bare g s . If all results are 
with bare g s , the present Hamiltonian gives the small- 
est r.m.s. deviation. We should also note, on the other 
hand, that the quenching factor obtained here has some 
ambiguity as the dependence of the r.m.s. deviation on 
the value of gf^/g s is quite modest. 



Table Q] presents the Gamow- Teller transition rates 
B(GT). The B(GT) values can be extracted from ex- 
perimental log ft values with the equation, 

ft = 6147 

1 (g A /g v fB(GTy W 

where 6147 is from Ref. [44[, g A and gy arc the axial- 
vector and vector coupling constants, respectively. For 
beta decays, we use bare gA/gv — —1-26 df|. The cal- 
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FIG. 15: (Color online) Similar to Fig. QU but for e p = 1.3 
and e„ = 0.5. 



culated results are with <7^ which is from ^-square fit- 
ting of these B{GT) values. The (g e / f /g A ) = 0.72,0.68 
and 0.64 for the present 2hu>, the present Ohuj and WBP, 
respectively. The (g^ /9a) value for WBP is very close 
to the commonly used value 0.60 [3l|. n Li and 15 C are 
weakly bound with 0.325 and 1.218 MeV neutron separa- 
tion energy, respectively. The protons in their daughter 
nuclei from j3 decay are well bound. Halo or skin effects 
need to be included which is not included in calculations 
with harmonic oscillator bases. The overlap between re- 
lated neutron and proton orbits is calculated in Woods- 
Saxon bases to modify the B(GT) of these two nuclei. 
All these three calculated results, WBP, the present in 
2hoj and OHuj, are modified by the halo or skin. More 
details can be found in Ref. Q. 

The present 2fuj results improve most of the B{GT) 
values compared with the present Ohu and WBP results. 
In order to show the difference between calculations 



■3 - 



c 

CD 

E 

02 
o 

H — ' 

CD 

D) 1 
CO 

E 





■ 


Expt. 


i — i — i — i - 

(a) 




□ 


WBP 




V 








A 


present 

A 




ft 




■ 


- 


ft I 


1 


m 
B 






■ 


A 

•!!,! 





20 O 19 O 18 O 17 O 16 O 15 O 13 O 19 N 18 N 17 N 16 N 15 N 14 N 13 N 12 N 



,3 - 



CD 

E 

o 2 
_o 

"-I—* 

CD 

£1 

CO 

E 



■ 


Expt. 




(b) 


□ 


WBP 






V 


present Oftco 


1 




A 


present . 


\ 


8 


* 






* 




I i 
' 8 


■ 

□ 


1 



17 C 15 C 14 C 13 C 11 C 17 B 15 B 14 B 13 B 12 B 11 B 10 B 



FIG. 16: (Color online) Magnetic moments calculated with 
WBP in Ohuj, present in both Ohio and 2ftw model spaces, 
compared with experimental data [34l ]. All magnetic moments 
are for the ground states except for 2^ in 20 O and ls O, 3j~ in 
16 and 14 C. 



and observed values, we present B(GT) exp / B{GT) ca i in 
Fig. HHJ B(GT) exp /B(GT) ca i from the present 2hu cal- 
culation is very close to unity except for 14 C and 14 0, 
and the second transitions in 15 C. 14 C and 14 are the 
same in the present isospin symmetric Hamiltonian. The 
abnormally long lifetime of 14 C has been a long-standing 
theoretical problem :46|. The present Qfvjj and WBP re- 
sults also fail in describing B{GT) of 14 C and 14 0. The 
reason is that two main components of the transition are 
almost all canceled in 14 C [4y|. It is hard to describe the 
cancellation exactly in interactions determined by consid- 
ering all nuclei nearby. In case of the second transition 
from 15 C to 15 N, the reason is similar, that is, three com- 
ponents are canceled resulting in a rather small value. 

Similar to the discussion in m.m., the r.m.s. deviation 
of calculated B(GT) values from the experimental ones 
is presented in Fig. It is clearly seen that both the 
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FIG. 17: (Color online) Root-mean-square deviation between 
calculated magnetic moments (m.m.) and the experimental 
ones, \m.m.(Exp.) — m.m.(Cal.)\, as the function of <?s /ff s - 



FIG. 19: (Color online) Root-mean-square of \B(GT) ca i 
B(GT) exp \ as the function of fl£ /<?A- 



VIII. SUMMARY 
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FIG. 18: (Color online) Ratio of the observed B(GT) values 
over the calculated B(GT), B(GT) exp / 'B{GT) cah in nuclei 
listed in Table UJ 



r.m.s. deviation and the quenching get smaller when the 
model space is enlarged. 



In the present work, we present a systematic study 
of boron, carbon, nitrogen and oxygen nuclei in full 
psd model-space with a newly constructed Hamiltonian. 
While some former Hamiltonians, such as PSDMK, WBP 
and WBT, are constructed in — lhui model space, we 
include 2 — 3Huj excitations in the present work. The 
present Hamiltonian is based on Vmu an d have four 
parts, (pp\V\pp) fromSFO, (sdsd\V\sdsd) from SDPF-M, 
(psd\V\psd) and (pp\V\sdsd) from Vmu plus spin-orbit 
force. We optimize the central part of Vmu while the 
tensor force in Vmu and the spin-orbit force are kept 
unchanged. The central force in (psd\V\psd) is 30% 
of Vmc/'s stronger than that in (pp\V\sdsd), while the 
strength of these two parts are the same in WBP. The 
SPE of the five orbits are also modified. More details of 
this Hamiltonian are explained in the text. 

The present Hamiltonian can reproduce well the 
ground-state energies, drip lines, energy levels, electric 
properties and spin properties of psc?-shell nuclei. Espe- 
cially, we can describe the drip lines of carbon and oxygen 
isotopes and spins of the ground states of 10 B and 18 N 
where WBP and WBT fail. The inclusion of 2fwj excita- 
tions is important in describing such properties because a 
part of mass-dependent effect in WBP and WBT is nat- 
urally included when we include 2ftw excitations. The 
effective operators become closer, in general, to bare op- 
erators when we enlarge the model space. We note that 
constant and smaller effective charges work quite well 
in the present study, which may attributed to the small 
size of the 4 He core also. The contribution coming from 
2Hu) excitations are investigated by comparison to Ohu 
calculations, suggesting that the present model space is 
still insufficient to reduce effective charges almost to zero. 
More systematic study is needed in a model space larger 
than psd and more hu excitations, especially for 4pAh 
excitations from p to sd shells. 
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It is also examined whether tensor force and spin-orbit 
force can be kept unchanged in full shell-model calcula- 
tions. Shell-model calculations without the modification 
of the strength of these two forces are found to be success- 
ful in the description of a wide range of psd-sheW nuclei. 
It is interesting to do more work on applying the present 
tensor and spin-orbit forces to shell-model calculations in 
other region of nuclei. 
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